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ABSTRACT 

An increasing body of data reveals a one-to-one linear correlation between galaxy halo mass and the 
total mass in its globular cluster (GC) population, Mqcs ~ M 103±0 03 valid over 5 orders of mag¬ 
nitude. In this paper we explore the nature of this correlation for galaxies of different morphological 
types, and for the subpopulations of metal-poor (blue) and metal-rich (red) GCs. For the subpopu¬ 
lations of different metallicity we find Mccsiblue) ~ M °- 96±003 anc j Mccsifed) ~ M^' 21±0 ' 03 with 
similar scatter. The numerical values of these exponents can be derived from the detailed behavior 
of the red and blue GC fractions with galaxy mass and provide a self-consistent set of relations. In 
addition, all morphological types (E, SO, S/Irr) follow the same relation, but with a second-order 
trend for spiral galaxies to have a slightly higher fraction of metal-rich GCs for a given mass. These 
results suggest that the amount of gas available for GC formation at high redshift was in nearly direct 
proportion to the dark-matter halo potential, in strong contrast to the markedly nonlinear behavior 
of total stellar mass versus halo mass. 

Of the few available theoretical treatments that d irect l y discuss the formation of GCs in a hierarchical 
merging framework, we find that the model of iKravtsov fc Gnedinl (120051) best matches these obser¬ 
vations. They find that the blue, metal-poor GCs formed in small halos at z > 3 and did so in nearly 
direct proportion to halo mass. Similar models addressing the formation rate of the red, metal-richer 
GCs in the same detail and continuing to lower redshift are still needed for a comprehensive picture. 
Keywords: dark matter — galaxies: formation galaxies: fundamental parameters — galaxies: halos 
— galaxies: star clusters: general — stars: formation 


1. INTRODUCTION 

A remarkably simple empirical correlation has emerged 
between two global properties of galaxies that belong 
to early phases of galaxy formation: these two quanti¬ 
ties are Mh (the total halo mass of a galaxy including 
all visible and dark matter in its potential well); and 
Macs (the total mass in its globular cluster system). 
Recent stu dies linking globular clusters (GCs) with halo 
mass (e.g. iBlakeslee et al.lll997t iSnitler fc Forbesl 120091: 

iGeorgiev et al.11201(1 iHudson et al.112014 among oth¬ 
ers; see Section 2 below for a more complete litera¬ 
ture survey) have shown that the simple ratio r/ = 
Mccs/Mh is virtually constant across 5 orders of mag¬ 
nitude in galaxy luminosity or mass. This behavior is 
quite unlike the strongly nonlinear shape of the more 
commonly k nown correlation of to t al stellar mass M+ 
with Mh Ce.g.lLeauthaud et, al.ll2012 : IHudson et al .11 201 Til 

iVelander et al.ll20l4 IKravtsov et al.ll2014ll . On strict ob¬ 

servational grounds, the new result indicates that there 
is at least one stellar subpopulation that formed in direct 
proportion to galaxy halo mass. 
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In our previous discussion (IHudson et all 120141 here¬ 
after HHH14 or Paper I), we constructed a new cal¬ 
ibration of = Mccs/Mh by using a recent cata¬ 
log of 419 GC system s (GCSs) in all types of galax¬ 
ies ((Harris et al.lhofi! hereafter HHA13). The GCS 
data were combined with a single, homogeneous calibra¬ 
tion of M* versus Mh founded solely on th e extensive 
CFHLens weak len sing program (Vela nder et al.l 12014 
IHudson et al.ll20T5l) . Measurements of M* in turn were 
determined from AT—band (or V —band in cases where 
K is unavailable) luminosities of the galaxies along with 
appropriate color-dependent mass-to-light ratios. The 
resulting correlation, built from 307 galaxies covering all 
galaxy types and all environments, confirms that Macs 
varies almost exactly one-to-one with Mh- The residual 
scatter is ±0.35 dex, and arguments taking into account 
the different weights of the individual Macs values and 
an estimated ;> 0.15 dex uncertainty in Mh suggest that 
the intrinsic scatter of the correlation may be no larger 
than ±0.2 dex. 

In HHH14 we treated all galaxies alike, and also treated 
their GCs as a monolithic population regardless of age or 
metallicity. In the present paper, we pursue the empirical 
correlation further by (a) subdividing the GC samples in 
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each galaxy into their astrophysically interesting “red” 
and “blue” metallicity groups, and (b) considering sub¬ 
groups of different morphological type (ellipticals, SO’s, 
and S/Irr types). We then look for any differences among 
these subsets that may provide additional clues for GC 
formation. 

In Section 2, we provide a more thorough overview 
of the observational and theoretical background for this 
topic, while in Section 3 we discuss some important fea¬ 
tures of the database. In Section 4 we describe the 
Macs ^ Mh correlation broken down by galaxy type and 
by GC metallicity groups. In Section 5 we discuss pos¬ 
sible interpretations, and in Section 6 we summarize the 
results. 

Readers interested primarily in the newer results and 
the analysis might find it useful to skim Section 2 on the 
background of the observations and relevant theory, skip 
Section 3 on the treatment of the database, and then go 
to Sections 4 and 5 for the results. 

2. BACKGROUND ON THE GCS -M H CORRELATION 

The idea that the GC population might scale directly 
with the total mass in its host galaxy has arisen several 
times over the past two decades. However, the previous 
literature on this topic is sparse enough that a fairly com¬ 
plete assessment of it can be made. To date, it has been 
an observationally driven subject. Thus we look first at 
the previous papers drawing on empirical data, and next 
at the (even fewer) papers directly addressing theoretical 
motivations for the correlation. 

2 . 1 . The Observational Material: Growth of the 
Database 

The first quanti t ative appearance of th e cor relation is 
inlBlakesled (|1997t l. lBlakeslee et al.1 (11997H . and lBlakesleel 
(119990 . who measured the GCSs in a sample of two 
dozen giant galaxies in Abell clusters. Several of these 
are BCGs (Brightest Cluster Galaxies), dominant ob¬ 
jects within their environments. Blakeslee and colleagues 
suggested in essence that the total GC population Ngc 
scaled linearly with the total (baryonic + dark) mass of 
their host galaxies as measured by halo X-ray tempera¬ 
ture or velocity dis persi on (see Fig. 15 and accompany¬ 
ing discussion of iBlakeslee et al.l 119970 . To quote from 
their paper, “These observations might be explained by 
a model in which the GCs ... formed early on and in 
proportion to the available mass”. However, their sam¬ 
ple was restricted to a small number of giant galaxies in 
unusually rich environments and at that stage the gen- 
e rality of the s tatem ent was not clear. 

iMcLaughlin (jl99 9) pursued a roughly similar line with 
a larger sample of 97 galaxies covering the size range 
from dwarfs to giants, focussing instead on the ratio of 
Nqc to the total estimated baryonic mass Mf, ary in a 
galaxy. In dwarf galaxies a large fraction of the initial gas 
mass has been lost due to feedback, and in supergiants, 
a large fraction of the initial baryonic mass is now in X- 
ray-luminous halo gas. If we add the extra assumption 
that the initial value of Mbary is proportional to Mh, 
i.e. Nqc ~ Mb ar y(init) ~ Mh, the argument connecting 
GCs to halo mass becomes similar in nature. Lastly, 
iKavelaarsi (|1999T ) used 11 large galaxies to show that Ngc 
scaled well with the dynamical mass of the galaxy within 


~ 20 kpc, as calculated from the velocity dispersion of 
the GCS itself. 

A decade later, a n ew t h read o f thi s empirical argu¬ 
ment was started bv lSoitler et al.l ( 2008 1 using a selected 
sampl e of 25 lar ge an d nearby galaxies. For this sample, 
as in IBlakeslee! (11997 1 they showed that the ratio of GC 
number to halo mass (Ngc/M h) was nearly constant, 
but now extended to g alaxies over a wider ran ge of envi¬ 
ronments and masses. iSpitler fc Forbesl (l2009ll discussed 
similar material for a larger sample of ~ 100 galaxies, and 
changed the thrust of the argument from GC numbers 
to GCS mass versus Mh . This difference is not trivial, 
because mean GC mass ( Mcc ) in creases sy s temat ically 
with galaxy luminosity; see IVillegas et al l (|2010H and 
HHA13. The calibration of Mh in Spitler & Forbes used 
a combination of two studies for weak-lensing masses of 
isolated galaxies, and (for lower-luminosity galaxies) con¬ 
ditional luminosity function analysis. 

iPeng et all (120081 1 used the GCS data from their HST 
survey of ~ 100 Virgo cluster galaxies to discuss both 
the luminosity fraction (denoted Sl) and mass fraction 
(denoted Sm) of GCSs per unit stellar mass M*. From 
these, they showed that the assumption Ngc ~ Mh can 
lead to the well known ve ry non-linear tr e nd of Sm ver¬ 
sus galaxy stellar mass. iGeorgiev et al.1 d 2010 H assem¬ 
bled Nqc totals for a still wider range of galaxies drawn 
from their own work on dwarf galaxies and from sev¬ 
eral other sources in the literature, along with galactic 
stellar and gas masses. They then defined the specific 
GCS mass as Sm = Mccs/Mbary where n ow Mhn.ru = 
M±+M gas , analogous to the ratio used by iMcLaughlinl 
(1999). The distin c tion b etween this definition and the 
one bv IPeng et al.l (12008H (who used Sm — Mgcs/M *) 
is significant for very gas-rich dwarfs and for giant ellip¬ 
ticals with massive amounts of X-ray halo gas. To trans¬ 
late Sm into 77 = Mccs/Mh and to calibrate Mh, they 
used a combination of dynamical mass measurements of 
dwarfs and galaxy group dynamics, though in a partially 
model-dependent way. For their sample of galaxies they 
find roughly Gaussian-shaped distribution functions for 
77 with small mean differences by galaxy type. 

HHA13 compiled all the Nqgs data available in the 
literature (419 galaxies at time of publication, contain¬ 
ing many more galaxies at high luminosities than in pre¬ 
vious discussions). This material was used to generate 
newly constructed values of Mgcs , explicitly accounting 
for the dependence of mean GC mass on galaxy lumi¬ 
nosity. They combined these with calculated dynamical 
masses for the galaxies, M ( j yn = const, a^R e . Since Md yn 
is computed at the effective radius R e within which stel¬ 
lar matter dominates, it is a reasonable proxy for Mbary 
or M*. Like the previous authors, they find that the mass 
ratio (Mgcs/M dyn) correlates nonlinearly with galaxy 
luminosity in a way that is consistent with a constant 77 
across the full range of galaxies. 

In HHH14 we took two additional steps: (a) the near- 
infrared luminosity of each g alaxy Mk was u sed along 
with mass-to-light ratios from lBell et al.l ( 2003 ) to repre¬ 
sent M*, and (b) M* was translated to Mh with a single, 
homo geneou s calibration based entirely on weak lensing 
data (|Hudson et al.ll2015fl. 

More recently. iFbrte et all (l2014f) have taken 67 galax¬ 
ies from the Virgo HST/ACS survey, divided them into 
9 luminosity bins, and constructed a ‘fiducial’ galaxy for 
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each bin that represents the average properties of the 
galaxies in the bin. They use these fiducials to derive the 
mean trend of Nqc per unit dark-matter halo mass for 
both the blue and red GC subsystems, finding nonlinear 
trends in opposite senses. Their calibration of Mh fol - 
lows the interpolation functions of lShankar et al. ( :200b ) 
based on occupation index and numerical simulations. 

A question of special interest is the presence of intra¬ 
cluster GCs within rich clusters of galaxies, which seem 
to belong to the overall potential well of the galaxy clus¬ 
ter as a whole and center on the BCG. T hese types 
of GC s have been foun d in Virgo ()Lee et all 1201 01 1 and 
Coma (|Peng et al. 11201 ll) and appear to be predominantly 
metal-poor, co nsistent with an origin from stripped satel¬ 
lite galaxies. iDurrell et al.1 (|2014f ) present new results 
from the NGVS (Next Generation Virgo Cluster Sur¬ 
vey) for the GC populations across the entire Virgo re¬ 
gion. They derive an estimate of 77 applicable to an entire 
galaxy cluster, where a large fraction of the total num¬ 
ber of GCs includes intracluster light. In this case, the 
dark-matter mass in Virgo is derived from X-ray halo gas 
temperatures and satellite dynamics. 

In Table 1, we summarize the different estimates of 77 
published in these observationally based papers and the 
number of galaxies used in each. The mutual agreement 
for 77 is within a factor of two. The source of the differ¬ 
ences from one study to the next is partly the steadily 
increasing size of the GCS databases used by each study 
in turn, but more importantly it is due to the differ¬ 
ences in the detailed assumptions about key parameters 
in the problem: the conversion of Nac into Mac Si the 
mean GC mass, the method for defining M*, and the 
conversion of M+ into Mh . These various assumptions 
are summarized in the last three columns of Table 1. 

2 . 2 . Theoretical Interpretations 

The history of theoretical models rel evant to the 
Macs — Mh correlation is even s horter. fPicman d et al.l 
( 2005) and IMoore et al.1 (12006H used high-resolution 
cosmological N-body simulations to conclude that the 
amount of material collected in high —a peaks at the very 
earliest stages of star formation (z ~ 12 ) should increase 
with the virial halo mass, roughly as M^ 2 (see Fig. 8 of 
Moore et al.). These peaks should be the places where 
proto-GCs would form. By hypothesis, then the oldest 
metal-poor (blue) globular cluster population should rise 
with halo mass, though perhaps nonlinearly. Their pic¬ 
ture, however, is closely linked to the role of reionization, 
which they and others b efore that (|Beaslev et alJl200% 
Santos 2003; Bckki 2005) suggested may have truncated 
GC formation at a slightly later stage. In addition, it is 
not clear how the metal-rich (red) GCs should fit in to 
this picture. 

Nearly simultaneously, (Kravtsov fc CnedTni (120051 . 
hereafter KG05) used a high-resolution hydrodynamic 
simulation to explore the formation of a typical Milky- 
Way-type galaxy, covering the time period from z = 11.8 
to 3.35 (set by limitations of computer time). This range 
covers the formation of the oldest (metal-poor, blue) 
GCs but unfortunately does not continue long enough 
to reach the [Fe/H] > — 1 range of the metal-rich, red 
GC population (see their Figure 8 ). The minimum cell 
size of their models is ~ 40 pc and thus individual GCs 
are not resolved. Instead, within each mini-halo of the 


merger tree, they identify the sites of GC formation as 
the densest cores above a threshold density of 1 M 0 
pc ~ 3 in the emerging GMCs (giant molec ular clouds) 
(cf. (Harris fc Pudritdll994 lKrulissenll2014 D. The input 
physics contains both internal SNe heating and exter¬ 
nal cosmic UV background. The time evolution of their 
model galaxies typically shows frequent bursts accompa¬ 
nying mergers and infall of gas-rich satellites, and rapid 
buildup of a massive central galaxy up to the end of the 
model run. The spatial distribution of the newly formed 
GCs ends up more centrally concentrated than the over¬ 
all halo potential, but most intriguingly for our purposes, 
the total mass in the GCs is found to track the masses 
of their host subhalos as Macs ~ M^' 13±0 ' 08 . This cor¬ 
relation is drawn from their simulated halos over a mass 
range 6 x 1O 9 M 0 <; Mh ^3x 1O 11 M 0 . 

In essence, both the Moore et al. and KG05 models 
should be compared more strictly with only the blue 
metal-poor GC population that emerged at higher red- 
shift. We emphasize this comparison in Section 4 below. 

In another cosmological simulation directe d towards 
finding GCs at formation, iBekki et~. (12008H identified 
proto-GCs as the central particles within virialized sub¬ 
halos. The simulation covers a large enough cosmological 
volume that many model galaxies can be inspected, thus 
building up GCS properties as functions of galaxy size. 
They find the peak rate of metal-poor GC formation to 
be at z ~ 8 — 6 and for metal-rich GCs at 2 ~ 4, in rea¬ 
sonable agreement with Milky Way GC ages (see above). 
For the ratio of Nac to halo mass, they find a shal¬ 
low, smooth increase of about a factor of two going from 
dwarf galaxies up to giants corresponding approximately 
to Nac ~ M 7 ) 3 , though with differences depending on 
galaxy morphology. A dependence this steep would not 
match the current data. 

Recentlv lKatz fe Ricottil (|2014T ) have modelled the pos¬ 
sible dependence of 77 on redshift (time of formation) by 
tracing the effects of an assumed 77 ( 2 ) on the present- 
day observed properties of the GCS. Given the amount 
of dynamical evolution and cluster destruction since for¬ 
mation, they deduce an increase of a factor of several at 
higher z, and also a hig her 77 in lo wer-ma ss halos. 

In both HHA13 and iKruiissenl (|2014f) . the empirical 
result 77 = const is viewed essentially as a product of two 
ratios, 77 ~ (M+/Mh) ■ (Ngc/M*). Each of these ratios 
varies nonlinearly with M*, but in opposite senses, so 
that when multiplied together their nonlinearities essen¬ 
tially cancel out. If so, GC formation might be largely 
disconnected from the dark matter halo and the unifor¬ 
mity of 77 is something of an accident. However, this view 
would give no explanation for the detailed way each ra¬ 
tio varies with galaxy mass, or why they should be so 
exactly complementary. 

An essential piece of the puzzle is the role of feed¬ 
back external to the proto-GCs during early star forma¬ 
tion. In dwarf galaxies, internal feedback including SNe, 
stellar UV emission, stellar winds, and external heat¬ 
ing from reionization may remove much of the original 
gas and truncate field-star formation. For giant galax¬ 
ies, feedback through AGNs and infall heating will simi¬ 
larly prevent much of the original gas in the DM poten¬ 
tial from forming stars. However, the extremely dense, 
small proto-GCs may have been much less affected by 
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Table 1 

Measurements of n = Macs /M^ 


Source 

n 

v Cio- 5 ) 

{Mqc) 

M* 

M h 

Blakeslee etal. 1997 

21 

17: 

2.4 x 10 s M q 

L v 

X-ray gas 

Spitler et al. 2008 

25 

3.2 

2.4 x 10 5 M© 

Lk,Lv 

weak lensing 

Spitler & Forbes 2009 

100: 

7.1 

4 x 1O 5 M 0 

Lk i Lv 

lensing + conditional LF 

Georgiev et al. 2010 

200: 

6 ± 1 

1.69 x 10 5 M 0 

L v + gas 

model, dynamics 

Harris et al. 2013 

255 

6.0 

function of My 
function of My 

alRe 

Lensing + abundance matching 

Hudson et al. 2014 

307 

3.9 ±0.9 

Lk,Lv 

weak lensing 

Durrell et al. 2014 

- 

2.9 ±0.5 

2.4 x 10 s M e 

L v 

X-ray, satellite dynamics 

This study 

175 

3.4 ±0.4 

function of My 

Lk,Lv 

weak lensing 


COLUMNS: (1) Literature source for the r]— value. 

(2) Number of galaxies used. 

(3) Derived value of mass ratio r). 

(4) Assumed mean GC mass needed to convert Nqq into Mqcs • For three studies (HHA13, HHH14, and the present work), both the 
mean GC mass {Mac) and the dispersion of the Gaussian GCLF are functions of galaxy luminosity My. Details of the calculation are in 
HHA13. 

(5) Method for deriving stellar mass M* ~ Mi lory . 

(6) Method for converting stellar mass to DM halo mass. 


any of these forms of feedback and more nearly reflect 
the original gas supply. As suggested before, the essential 
two-step link from halo mass to GCS mass is therefore 
Mqcs ~ Mbary{init) ~ Mh- As also seen in the theoret¬ 
ical simulation studies listed above, the role of the DM 
halo is to set up the conditions where dense gas collects 
together and the proto-GCs can form. 

As noted above, several previous discussions of GC for¬ 
mation are intimately connected with the role of cosmic 
reionization. In most of these papers, reionization is in¬ 
voked to truncate the formation of the earliest, metal- 
poor GCs and thus help produce th e observed b imodal 
metallicity distribution of GCs (e.g. iBeaslev et al.ll2 00% 


metaincity distribu tion or bus (e.g. m casicv ct ai. zuuz; 

Santosll2003tIBe kki 2005; Spi tler et al.ll2012tlMoran et al.1 


20141 ). This approach would require that most of the 
metal-poor GCs should have formed before the reioniza¬ 
tion peak, i.e. in the regime z )> 10. A second view 
is that reionization triggered GC formatio n by shocking 
protocluster clouds into collapse (lCenll200H ) . Still a third 
view is that the reionization epoch is largely irrelevant 
to GC formation, either because the proto cl uster c l ouds 
and their GMCs are self-shielded (IMuratov fe Gnedinl 
120101 : iBarkana fe Loebll200H iShapiro et al.ll2004D , or be¬ 
cause mo st GCs formed after the main epoc h of reioniza¬ 
tion ie.g. IBekki et alJfeOOSt lKruiissenll20l4 and KG05). 
KG05 suggest particularly that the conditions needed for 
GC formation (namely, growth of very dense cores em¬ 
bedded within massive, cooled GMCs) could only accu¬ 
mulate in a large-scale way after reionization and should 

have peaked around z ~ 5. __ 

Reionization peaked at z ~ 1 1 (jHinshaw et all 120131 
iPlanck C ollaboration et al.ll2014f) and was c omplete near 
z ~ 6 ( Fan et al.l 120061 : iZahn et aDl2012t |j ensen et all 
I2014D . By comparison, direct measurements of GC ages 
that do not depend on cosmological models find ages 
for Milky Way GCs distributed rather evenly from 11 
to 13 Gy , equivalent to a redshift range z — 2 5 — 7 
(see, e.g. iMarm-Franch et all 120091: lHansen et all [20131 


iLeaman et all 120131 iVandenBerg et al.l 120131 for recent 


analyses). With the standard techniques of isochrone fit¬ 
ting to deep color-magnitude diagrams for GCs, internal 
precisions of ±0.5 Gy are now achievable, though uncer¬ 
tainties of ~ ±1 Gy in the absolute ages remain due to 
small uncertainties in the different stellar model calcu¬ 
lations, the adopted GC distance scale, the adopted GC 
metallicities, a nd details of isochrone fi tting procedures 
(see especially IVandenBerg et all 120131 for a thorough 
discussion). 

If the GC age measurements summarized above are 
correct, the oldest GCs formed clearly before the peak 
epoch of star formation in the universe at 1 z <, 3 
(iBehroozi et al.lf2T)13L iMadau fc Dickinsonl[2014l ). but af¬ 
ter reionization had largely finished (though recognizing 
the caveat in the absolute GC age calibrations). The rea¬ 
son for bringing in reionization was to produce the uni¬ 
versal GC bimodality, but this may also be unnecessary. 
An alternate and perhaps more attractive explanation for 
bimodality emerges from the fact that the merger tree for 
a large galaxy contains very few intermediate-sized halos 
in which intermediate-metallicity GCs would form; most 
mergers are strongly unequal-mass o nes w here dwarfs 
are accreted by the central galaxy dCote et al.l 119981 : 
IMuratov fe Gnedi32010HToninill2013l) . These arguments 
combined favor the third option that cosmic reionization 
and other forms of feedback had far less influence on GC 
formation than on the lower-density field star population. 

Though it is at early stages, the theory available so far 
suggests we should look more closely at the Mqcs — Mh 
correlations for the easily identifiable GC subpopulations 
(blue, red) and also host galaxy type. We do this in 
section 4 below, after some comments about the nature 
of the GCS database itself. 

3. THE GCS DATABASE AND TESTS OF SYSTEMATICS 

In the present paper we continue to work with the data 
in the HHA13 catalog of globular cluster systems. Be¬ 
cause the catalog is assembled from more than 100 in¬ 
dividual GCS studies in the literature, inhomogeneity is 
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Figure 1. Correlation of the total mass in globular clusters, 
Mgcs 5 versus the total halo mass Mh of a galaxy. Masses are 
in units of Mq. The errorbars on each point show the uncer¬ 
tainties due to the measured total GC population in each galaxy. 
Magenta open circles are galaxies for which the GC population was 
measured through the SBF technique (see text). The solid line is 
the least-squares solution for the “All” sample listed in Table 2 
and has a slope of unity. The shaded region at bottom shows the 
range of low-mass galaxies that would on average be expected to 
have less than 1 GC (and thus would not appear in our database). 
At lower left, the gr e en lin e shows the theoretical prediction from 
I Kravtsov 8z Gnedinl H20051) over the halo mass range covered by 
their model. 

undoubtedly present and may be responsible for some 
of the observed scatter in the basic correlation of Mac s 
versus Mh (e.g. iZaritskv et al.ll20Tl : lUsher et al.ll2018fl . 
However, the estimated Ngc values in most cases are de¬ 
rived from a simple and highly standardized method now 
used in almost all GCS studies. In brief, this “standard 
method” is: 

• Begin with imaging of the target galaxy that is as 
deep as possible, as wide-held as possible, and uses 
as many filters as possible, 

• Carry out photometry across the imaging held 
down to some objectively defined limiting magni¬ 
tude (usually the 50% completeness level as deter¬ 
mined by artificial-star experiments). 

• Subtract the background level (held contamina¬ 
tion). The GCS shows up as a statistical excess 
of starlike or near-starlike objects, centered on the 
galaxy and showing up above the uniform back¬ 
ground of ‘held’ objects. The held contamination 
can be reduced by a number of methods including 
elimination of extremely red or blue objects, ob¬ 
jects that have structures marking them as faint 
background galaxies, and objects that have very 
different scale sizes than GCs at the distance of the 
target galaxy. For HST-type resolution, the GCs 
may be partially resolved and thus not starlike. 


Figure 2. Correlation of the total mass in globular clusters, 
Mgcs 5 versus the total halo mass Mh of a g alaxy, now for 8 4 
galaxies in the HST-based Virgo Cluster Survey )|Peng et all20081 ). 
The dashed line shows the best-fit solution for all galaxies, from 
Figure 1, while the solid line shows the solution from only the 
Virgo members. The green line at lower left is the KG05 model 
line as in Figure 1. 

• Calculate the residual total measured GC popula¬ 
tion in the field after background subtraction. Cor¬ 
rect this total for the unobserved outer-halo portion 
beyond the field of view of the study, normally us¬ 
ing the power-law or Sersic-law radial distribution 
of the GCs as measured directly from the data. 

• Finally, correct the total to include GCs fainter 
than the photometric limit of the data. Here, 
we gain a large advantage from the fact that the 
GC luminosity distribution (GCLF) has a stan¬ 
dard lognormal shape in galaxies of all types and 
sizes, a long-established result resti ng o n extreme ly 
strong observational footing dVillegas et all 120101 : 
l.lorda.n et a.l.ll200'7aL lHarris et al.ll2f)14blb 

Excellent recent examples of the application of this 
standa r dized metho d can be found in IZaritskv et al. 
( 2014T): iBrodie et aLl (I2014T) : lYoung et all (12012ft : iHarri* 
(2009al): iCho et all (120121) to mention only a few. An 
extreme example is the uncovering of the GCS in the 
nearby giant elliptical NGC 51 28 holdi ng a system of 
~ 1000 GCs ([Harris et, al.ll2004aHbl . I20l3l . This target in 
some sense represents the most difficult combination of 
all the selection techniques listed above: the surround¬ 
ing field has large amounts of contamination from both 
foreground stars and faint background galaxies; spatial 
resolution is not as effective since the GCs are marginally 
resolved under good ground-based seeing; and the galaxy 
is so nearby that its GCS is spread across a large pro¬ 
jected area of the sky. Even so, reliable results for the GC 
total Ngc and GCLF have been ob tained, with several 
studies showing mutual consistency (lHarrisll201ClD . 
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The most prevalent sources for any inhomogeneity in 
the current catalog (HHA13) are likely to be field size 
and radial coverage, which differ qui te a bit from one 
individual study t o another (see also lllsher et al.1 1 201. 'll 
IBrodie et al .1120 l ill . The galaxies most affected by lim¬ 
ited field coverage are giant, nearby ones that have been 
observed with cameras with relatively small fields of 
view, such as the HST WFPC2, ACS, or WFC3 detec¬ 
tors. In some studies, careful work has b een done to as¬ 
sess the field corrections (e.g. lPeng et al.ll2008lh even cal¬ 
ibrated against wide-held ground-based imaging, but the 
resulting uncertainties in Ngc are inevitably larger than 
for studies drawing directly on wide-held imaging. For 
the limiting-magnitude correction, most modern studies 
reach or exceed the GCLF peak or ‘turnover’ point above 
which half the GC population is found, and the univer¬ 
sally consistent shape of the GCLF makes the resulting 
normalization to Ngc easy to carry out. 

The issue for the present paper in practice boils down 
to two questions: (a) does inhomogeneity systematically 
affect either the slope or zeropoint of the correlation; 
and (b) does it generate significantly increased scat¬ 
ter 1 ! An empirical test of both questions can be made 
with the study in the literature that probably has the 
strongest claim to internal homo geneity, the H ST-based 
Virgo Cluster Survey IVCS: IPeng et alll2008lh Here all 
target galaxies are at virtually the same distance and 
were observed with the same camera (ACS/WFC), the 
same exposure times, and the same filters. Even here, 
however, field size is still an issue because a smaller frac¬ 
tion of the total GCS falls within the field of view for 
larger galaxies. 

The basic correlation of Mgcs versus M /,. including 
all galaxies with measurements of both Ngc and Mf, is 
shown in Figure [TJ and for comparison Figure [5] shows the 
correlation for only the 84 Virgo members with both halo 
masses fr om HHH14 and Ngc numbers adopted from 
iPeng et all (1200811 PI Encouragingly, the best-fit solution 
from the Virgo data alone is not significantly different 
from the solution from all catalog data. However, the 
scatter is certainly lower (formally 0.30 dex for Virgo, 
versus 0.35 dex for all data). 

The entire catalog has many more galaxies than does 
the Virgo sample in the mass range Mh > 10 12 Mq, in¬ 
cluding all but one of the BCGs, and it is in this high- 
mass range that the scatter is perhaps most problem¬ 
atic. Though the inhomogeneity of the data for the larger 
galaxies may well be the cause of the increased scatter, 
it is hard to state this definitively since the Virgo cluster 
contains only a few such galaxies, so it may also be pos¬ 
sible that some of the scatter is intrinsic. A single, more 
homogeneous observational survey targeted specifically 
to a large sample of high- mass gala xies would provide a 
key test of this point (see lZaritskv et aL| [20144 as well as 
being valuable in its own right. 

For each distinct subset of the data we carry out a 
separate best-fit solution from a weighted fit of all data 
points. Specifically, we assume that data pairs ( Xi,yi) 

2 We note here that Virgo contains two supergiant ellipticals 
of similar visual luminosity, M49 and M87. M87 is more likely 
to be the Virgo BCG because it is sitting near the center of the 
large-scale Virgo potential well and is surrounded by an extended 
halo of hot X-ray ga s, intracluster field stars, and intracluster GCs 
(IDurrell et al.|[20141) . 


with measurement uncertainties (< T x ,ii&y,i ) are related by 

y = a + /3(x - (x)) (1) 

and then minimize the y 2 estimator dNovak et all 120061 : 
lHarris et al.ll2014at 

2 = V- {Vi -a-P(xj - (x))) 2 

X ^Ki + eD+Ptej + el)' U 

In this case, x = log Mh and y = log Macs- In Eq. [2j e x 
and e y denote any intrinsic or “cosmic” scatter present 
in the data and additional to the direct measurement 
uncertainties. Without extra knowledge of some kind 
about the data, we would not be ab le to gau ge both e x 
and e y independently (see lHarris et al.ir2014al) . 

Table 2 lists the results for all the datasets discussed 
in this paper. Following HHH14 we set a x = 0.15 dex for 
all points, whereas a y is determined by the uncertainties 
in Ngc as given in the catalog. For the initial solutions 
listed in Table 2 we set e x = e y = 0, lacking extra infor¬ 
mation that would specify otherwise. The last column 
of the Table gives the total rms scatter of the residuals 
in y. For all the solutions listed, the fits were performed 
iteratively removing any objects with y —residuals larger 
than 1.0 dex. The n —values listed in column 3 give the 
number of galaxies remaining for the final fit. 

The HHA13 catalog contains about 40 galaxies for 
which Ngc was measured with the SBF (surface 
brightness fluctuation) techniqu e ((Blakcslcc et alJHf)97t 
iMarfn-Franch fc Anariciol 120021 1 where the GCS is de¬ 
tected through the imprint of the unresolved GCs on 
the power spectrum of the integrated light, rather than 
the standard method outlined above. These objects are 
plotted separately in Fig. 1, and in Table 2 we list the 
solutions both including and excluding them. These two 
solutions are the same within the uncertainties of the co¬ 
efficients ( a , /?), but the SBF technique does not fall into 
the ‘standard method’ category and in what follows we 
exclude any data from this method. 

To investigate two other possible systematic effects, in 
Figure [3] we show the residuals from the Mgcs — Mh 
correlation plotted versus galaxy distance d , and versus 
galaxy luminosity My. The data shown are for objects 
with the most well determined GCS totals (< 0.15 dex 
in Ngc ) and also exclude objects measured only with 
the SBF technique mentioned above. The halo masses 
are determined from the integrated near-IR luminosities 
Ktoti transformed into M* and then into Mh through 
the calibration curve described in HHH14. Thus if the 
luminosities of large, nearby galaxies are underestimated 
because of aperture-size limitations, these galaxies would 
show positive residuals since their predicted halo masses 
would be too small. No obvious effect of this type shows 
up in Figure 3a. Similarly, in Figure 3b there is no 
strong trend wit h gala xy luminosity. As discussed in 
iSnitler fe Forbes! (I2009D and HHH14, for many of the 
BCGs particularly in large Abell clusters, Mgcs is likely 
to be underestimated because their GC populations are 
so spatially extended that they encounter field-of-view 
limits in most observational studies. Individual correc¬ 
tions for field size and intracluster populations must be 
made and are difficult to predict. To date, these specific 
corrections h ave been m ade on ly for the BCGs in Virgo 
(M87; see IDurrell et al.l 120141 ) and Coma (NGC 4874; 
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Table 2 

Observational Correlation of Macs with 


Galaxy Type 

GC Subset 

n 

«(±) 

/3(±) 

(x) 

a (rms) 

All Types 

All 

293 

7.706 (0.020) 

1.000 (0.020) 

12.3 

0.349 


Blue 

167 

7.405 (0.025) 

0.960 (0.027) 

12.2 

0.323 


Red 

145 

7.157 (0.030) 

1.209 (0.030) 

12.2 

0.355 


No SBF 

269 

7.598 (0.021) 

1.004 (0.019) 

12.2 

0.347 


f rP R known (“Best”) 

166 

7.564 (0.024) 

1.026 (0.025) 

12.1 

0.306 


Virgo only 

82 

7.118 (0.033) 

1.035 (0.043) 

11.7 

0.302 

E 

All 

171 

8.041 (0.027) 

0.960 (0.023) 

12.6 

0.348 


Blue 

97 

7.549 (0.035) 

0.945 (0.036) 

12.2 

0.340 


Red 

86 

7.321 (0.041) 

1.186 (0.041) 

12.3 

0.376 

SO 

All 

77 

7.627 (0.037) 

1.150 (0.049) 

12.2 

0.326 


Blue 

47 

7.253 (0.037) 

0.980 (0.056) 

12.0 

0.254 


Red 

38 

7.063 (0.045) 

1.304 (0.037) 

12.2 

0.280 

S/Irr 

All 

46 

7.086 (0.052) 

0.985 (0.078) 

11.8 

0.350 


Blue 

24 

7.171 (0.068) 

0.825 (0.085) 

12.0 

0.334 


Red 

21 

6.913 (0.070) 

1.107 (0.114) 

12.0 

0.322 


see lPeng et alJl20llh . Contrarily, Mh is not as affected 
because the prescriptions for defin ing halo masse s in cen¬ 
tral galaxies such as BCGs dHudson et al.ll20T5l) include 
their extended halos. The net result would be to leave 
negative residuals, which is what we see for many of the 
highest-mass galaxies in Fig. [D Most of these are miss¬ 
ing in the more highly selected sample of Fig. [3] In the 
only two cases (Virgo, Coma) where it has been possible 
to perform a wider-scale census of the GCs in the entire 
cluster, this negative bias is also reduced. 

In summary, we have no evidence that use of the entire 
GCS catalog is preventing us from finding the necessary 
global trends. With the BCG-related problem in mind, 
it appears worthwhile to use the entire GCS catalog of 
data to investigate other astrophysically interesting cor¬ 
relations, with GC metallicity and galaxy type. These 
are discussed in the next sections. 

4. THE MASS RATIOS VERSUS METALLICITY AND 
MORPHOLOGY 

4.1. Cluster Metallicity: Blue vs. Red 

In most galaxies, the GC metallicity or color distri¬ 
bution has a bimodal shape with the blue (metal-poor) 
and red (metal-rich) modes often partially overlapping. 
The population fractions are f re d = N(red)/N(tot) and 
fbiue = 1 — fred- With the assumption that the mass 
distribution functions for both blue and red GCs are the 
same, then we also have f r ed = Mgcs(red)/Macs (see 
Villegas et al.l l2010t lAshman et al.lll995t iTamura et al.l 
20061 : lLarsen et al.l2001[) . The blue mode is centered near 
(Fe/H) ~ —1.5 and the red mode near (Fe/H) ~ —0.5, 
though the peaks of both mode s show second-order in¬ 
creases with galaxy luminosity (iBrodie fc Strader! [20061 
iPeng et al.ll2006l : [Forte et al.ll2009D . The dividing line be¬ 

tween the blue and red modes is conventionally defined 
as the central minimum in the color or metallicity his- 



0 1 2 
log (d/Mpc) 



-16 -18 -20 -22 -24 

M V T (galaxy) 


Figure 3. Upper panel: Residuals (measured minus predicted) 
from the best-fit correlation of Mqcs versus M h, plotted versus 
galaxy distance. The sample shown includes only galaxies with 
well determined GC populations (uncertainties < 0.15 dex) and 
excludes any measured with the SBF technique (see text). Lower 
panel: Residuals plotted versus absolute visual magnitude of the 
galaxy. 


togram, which usually lies near [Fe/H] ~ —1.0 as it does 
in the Milky Way. 

We have inspected the literature sources for all galax¬ 
ies in the HHA13 catalog and wherever possible extracted 
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Figure 5. Mass ratio of the GCS to the halo mass rj versus halo 
mass, for the 175 selected galaxies with best data (see text). BCGs 
are shown as magenta open circles. The green line marks a running 
median value of ?/ calculated in bins of 12 galaxies each. 


Figure 4. Correlation of Macs with AT/, for the “best” subsample 
of 175 galaxies as defined in the text. BCGs are shown as the 
magenta open symbols at upper right, and the best-fit line has a 
slope of 1.03 ± 0.03. The green line and shaded region at bottom 
are as defined in Fig. [T] 


the fractions f re d ■ These values will be listed in the next 
catalog release. In many cases the individual studies al¬ 
ready list fred , usually derived from a double-Gaussian fit 
to the color or metallicity distribution. In a few papers, 
we have reconstructed them from the available data. For 
almost half the list, however, the red/blue ratios could 
not be determined because the photometric studies were 
carried out in single filters, or the sample size was too 
small, or the color indices were not precise enough to 
identify bimodality clearly. Any GCSs measured through 
only the SBF method are also excluded. The net result 
is a list of f re d values for 219 catalog galaxies, of which 
175 have calculated halo masses M /, and 14 are BCGs. 
We label this the “best” subsample of catalog data0 

The Macs ~ Mh correlation for this selected dataset 
is shown in Figure [U for which we find a best-fit solution 
Mqcs ~ M^ 03±0 - 03 . The residual scatter is 0.3 dex, sim¬ 
ilar to the quality of fit for the Virgo subsample (Fig. [2j 
but with the advantage that the higher-luminosity range 
of galaxies is more thoroughly covered. Interestingly, 
many of the BCG and other giant galaxies that had 
negative residuals in Fig. [Tj have dropped out of Fig. 0] 
Some of these large, distant galaxies had lower-quality 
data that prevented blue and red fractions from being 
determined reliably, and some were measured with the 
SBF technique that automatically excludes them from 
the red/blue plot. The distribution of 77 for these 175 
galaxies is shown in Figure [SJ The weighted-mean value 
of the mass ratio is (log??) = —4.47 ± 0.07 (listed in the 

3 The “best” sample is defined as the galaxies for which the 
source photometry was good enough to determine the red and blue 
GC fractions; though these are not necessarily the ones with the 
smallest nominal errorbars on Nqq. 


last line of Table 1). 

In Fig. [5l we also show a running median value of r) 
in bins of 12 galaxies each, where the successive bins are 
stepped through the list sorted by M h . The median does 
not deviate far from 77 ~ const, and we conclude for the 
present that the simple linear form Macs ~ Mh cannot 
be ruled out by the data. The range of most concern is 
probably the high-mass end Mh )> 10 13 Mq, where most 
of the points fall below the global average. However, too 
much significance should not be put on these biggest, 
highest-mass points for at least two reasons: first, a slight 
change in the adopted power-law slo pe of the conversio n 
between M* and Mh at the high end (IHudson et al 120151) 
would remove this mean offset entirely. Second, as men¬ 
tioned above and discussed by HHH14, the estimated 
GC populations for many of these tend to be underesti¬ 
mated because of the huge spatial extent of their GC sys¬ 
tems, which smoothly merge into wh atever intracluster 
GC population may be present fe.g. lDurrell et al.l [2014). 

The correlation of f re d with galaxy luminosity My is 
shown in Figure [ 6 j More luminous gal axies te n d to have 
higher fractions of metal-rich GCs dPeng et al.ll2006D . but 
with considerable scatter. Many dwarfs have only blue 
GCs present ( f re d = 0, at bottom of figure). The mean 
trend is very roughly given by f re d — 0.25 — 0.052(My + 
19), with residual scatter equal to ±0.13 rmsQ 

The large scatter is partly intrinsic, but inspection 
of the source literature also indicates that part of it 
arises from differing and incomplete field coverage be¬ 
tween studies, as discussed above. Because the bluer 
GCs are spatially more extended around the host galaxy 
than the red ones, any study covering only the inner re¬ 
gions of the halo will overestimate f re d unless the radial 
distributions are explic itly corrected for (these correc- 
tions are done, e.g., in iPeng et all 120081 : iHarrisl [200911 

4 If we use log M* on the x-axis instead of My, the distribution 
is very much the same and with the same scatter, but with fewer 
datapoints because the Ktot IR luminosities are not available for 
many galaxies. 
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Figure 6. Fraction of the GC population that is metal-rich (red), 
plotted versus host galaxy luminosity My. BCGs are plotted as 
open circles. The mean correlation given in the text is shown with 
the dashed line. 



log (M holo /M e ) log (M halo /M e ) 


Figure 7. Left panel: Total GC mass in the metal-rich GCs versus 
halo mass. Right panel: Total GC mass in the metal-poor GCs. 
The KG05 model prediction is shown as the green line at lower left. 
In both panels the best-fit solution from Table 2 is shown as the 
solid line, while the solution for all GCs combined is the dashed 
line. 

iDurrell et al.ll20fl) . Nevertheless, it is worth noting that 
in 86% of all galaxies in the sample, the majority of the 
GCs even from the uncorrected totals are the blue, more 



Figure 8. Correlation of the total mass in globular clusters, 
Mccs ■ versus the stellar mass M* of a galaxy. Symbol types 
are as in Fig. [7] The dashed line in each panel shows the expected 
relation for galaxies with a GC specific frequency Sc = 1 (see 
text). 


metal-poor ones, and even in many of the supergiants. 
Galaxies where the metal-richer GCs dominate are quite 
rare, and these few are not necessarily the most mas¬ 
sive galaxies. For the central giant galaxies particularly, 
many of their outer-halo GCs should be from captured 
dwarf satellites that mostly contained metal-poor GCs. 
The intraclustcr G C population considered by itself is 
largely metal-poor dLee et al .11 201(1 IDurrell et al.ll2014D . 

Figure [7] shows the results for Mccsired) = freduces 
and Mccs{blue) = fbiueMocs plotted as before versus 
Mh- For comparison, Figure [8] shows the same data plot¬ 
ted versus M*. In Fig. 0 note that the KG05 model 
line is displayed only for the blue-GC subgroup, since 
the metallicity range it covers does not apply to the 
red subgroup. Linear least-squares solutions for each 
are shown as the solid lines and listed in Table 2. For 
these solutions, we adopt an extra dispersion e y = 0.1 
given the uncertainties mentioned above in defining f re d- 
For M( 3 c>s(red) we deliberately exclude any dwarfs with 
fred = 0, to concentrate on the trends for galaxies mas¬ 
sive enough to have generated metal-rich GCs. Ear¬ 
lier versi ons o f the se trends using smaller sample were 
shown bvlSoitler et al.1 (120081) versus both Mh and Mb ar y, 
and bv lForte et all (120141 ) for Nqc versus Mh- Forte et 
al. used the Virgo galaxies binned by luminosity; we pre¬ 
fer to show all galaxies individually to gain a clearer idea 
of the scatter and relative quality of the correlation. 

As expected given the dominance of blue GCs in most 
galaxies (Fig. 02), the blue slope Macs(blue ) ~ M° 96 is 
only slightly shallower than the Macs ~ M^' 03 slope for 
all GCs. The steeper slope Mccs^ed) ~ Ml' 21 tracks 
the more rapid rise of the metal-richer subpopulation 
during the enrichment history of the galaxies concerned, 
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with a hint of nonlinearity as well in the mid-range of 
galaxy masses. Notably, the scatters for both subpop¬ 
ulations are similar to each other and to the scatter we 
found for the combined population (rms ±0.32 — 0.36 
dex). 

The galaxy mass at which f re d or Macs (red) rather 
suddenly plunges to zero is near Mh ~ 10 11 Mq or equiv¬ 
alently M* ~ 10 9 Mq. In physical terms, this thresh¬ 
old represents the minimum halo mass capable of reliably 
generating , and holding, GCs with metallicities [Fe/H] 
> - 1 . 

In Fig. [8] the dashed line shows the expected re¬ 
lation f or ga l axies with speci fic frequency Sn = 1 
(iHarris fe van den Berghl I1981H . where Sn = Nqc ■ 
10 “°'4( M v+i5) anc t Nqc has been converted to Mqqs 
following the relations in HHA13. Neither relation can 
be adequately matched by a single power-law slope over 
the entire range of galaxy size (this is the classic “specific 
frequency problem”). The red GCs are, however, closer 
to varying linearly with M* than are the blue GCs, as 
expected if their formation epoch was closer to the uni¬ 
versal peak of star formation around z ~ 2 that produced 
most of the galaxies’ stellar mass. There is a restricted 
intermediate-luminosity range near 10 logM* 10.7 
(blue GCs) or 9 logM* < 10.7 (red GCs) where it is 
possible to claim that GCS mass is directly proportional 
to stellar mass. This is precisely the range over which 
Sn is constant for all types of galaxies (see Fig. 10 of 
HHA13). As noted earlier, M* is not a tr ue prox y for 
total baryonic mass in ma ny galaxies (see [McLaughlin 
119991 iGeorgiev et al.ll2010l f and it would be desirable to 
assemble a more complete contemporary listing of Mb ary 
for the galaxies in the catalog. 

In short, we find that galaxy halo mass is a much sim¬ 
pler predictor of GC population than galaxy stellar mass, 
regardless of GC metallicity. 

The Kravtsov/Gnedin (KG05) model prediction is 
plotted in Figs. Q] and [7b at lower left. The agreement 
in both slope and zeropoint with the real data from the 
dwarf galaxies is striking. Most importantly, the model 
slope of 1.13 ± 0.08 is encouragingly close to the observa¬ 
tional result of 1.03 ±0.03 for the “Best” sample defined 
above. The zeropoint agreement should, however, be re¬ 
garded as more nearly fortuitous: the model data are the 
estimated protocluster (dense core) masses at formation, 
whereas the observations give the present-day values of 
Macs, which are likely to be ~ 2 — 3 times lower than 
their initial values after early mass loss and a Hubble 
time of dynamical evolution. In addition, the KG05 halo 
masses are evaluated at z > 3, while the observationally 
based numbers we use here for Mh are the halo masses 
at z ~ 0. 


4.2. Galaxy Type 

An entirely new correlation that can be derived from 
our database is the trend of Macs with Mh, plotted 
separately for three basic morphological types: ellipticals 
(—6 < T < —3), SO’s (—2 < T < 0), and spirals and 
irregulars together (T > 0). The data are displayed in 
Figure [9] and best-fit solutions listed in Table 2 as before. 
Here also the BCGs are not marked separately, though all 
of these belong to the E group. The scatter is smallest for 
the SO’s, likely because many of those are drawn from the 


1 1 1 1 ■ 1 1 ■ 1 1 

/ 

1 1 1 1 1 1 1 1 1 1 

... 

t*' 


/ 

E * • */» 

SO ' 

S/Irr ' 


?• 

/ 

•4-* 

• . /• * 

• • 

1 

* / 

4 

I 

:> ! > 

T 

V* 

• / 

/• 

r 

i 

_ •/ • _ 

•/. 

J* 

s*f 

/ 

£ 

f. 

• <• • • 

1 

7* 

1 

_1_ 

i 

i 

• • 

■ * 

.t- 

/. •* 

/ * 

/ 

/ 

— ™ • — 

/. 

- / * 

/ 

- / 

/ 

- r 

/ 

/ 

7 • 

/ 

7 

/ 

7 

. i i i , i i ,< i ,<- 

< i i i i i , i i i i >. 

• ImiIm.Iih 

10 12 14 

10 12 14 

10 12 14 


lo g ( M halo/ M s) 


Figure 9. Correlation of the total mass in globular clusters versus 
galaxy halo mass, for different galaxy morphological types. There 
are 179 data points for ellipticals (left), 78 for SO’s (center), and 
46 for spirals and irregulars (right). In each panel the dashed line 
shows the mean correlation for all types combined, from Figure 1. 


internally more homogeneous Virgo and Fornax surveys. 

In Figure [TO] we show the best-fit mean lines for all 
three galaxy types on the same graph, along with the 
solution for all galaxies combined for intercomparison. 
Each line is plotted over the range of masses occupied by 
the available data for that subgroup. Though the slope 
for the SO’s is highest by a noticeable margin (/3 = 1.15 ± 
0.05 versus 1.03 ± 0.03 for all galaxies), the difference 
is not strong and it is difficult to know what, if any, 
significance to ascribe to it. 

A potentially more important difference - though sim¬ 
ilarly small is that the mean lines for the S/Irr and E 
types are nicely parallel but offset by (0.18±0.06) dex. If 
taken at face value, this offset would mean that the late- 
type galaxies were roughly 60% as efficient as the early 
types in forming GCs. Again, however, this difference is 
just at the margin of detectability with the current data. 

An alternate explanation for the offset between the 
E/S0 and S/Irr mean efficiencies might be environmen¬ 
tal, in the sense that early-type galaxies are found more 
often (though not exclusively so) in richer groups and 
clusters of galaxies. Denser environments might logi¬ 
cally promote hi gher GC and star formation efficien cy at 
high z (see, e.g. iPeng et al.l 120081 : lKruiissenll2014fl . We 
will explore the environmental density parameter more 
quantitatively in a subsequent paper. 

5. DISCUSSION 

There is no guarantee that a new set of data will help 
decipher the history of the systems we are studying. But 
the remarkably tight and structurally simple link be¬ 
tween GCS mass and galaxy halo mass suggests to us 
that it can add to our understanding especially of GC 
formation and more generally their environments at time 
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Figure 10. Mean relations from Table 2 for three different galaxy 
morphological subgroups (E in red, SO in green, S/Irr in blue). 
The dashed line is the solution for all galaxies combined, as in the 
previous Figure. 

of formation. 

5.1. Galaxy Morphology 

The message from subdivision by galaxy type (Fig. 0 
seems straightforward. All major types follow much the 
same relation as defined by the entire sample, though 
with a possible second-order offset between ellipticals 
and spirals. The standard Hub ble se q uenc e was not 
well established till z <j 2 (e.g. IMortlock et al.1120131 
iFiacconi et al.l2015l) , clearly after the main era of GC for¬ 
mation. The simplest conclusion is that the majority of 
the GCs formed at an early enough time that their condi¬ 
tions of formation were not strongly sensitive to the later 
emergent galaxy type as we now see it. This evidence 
reinforces other striking similarities in GC populations 
that have long been known. Across all types and sizes of 
galaxies, GCs have (a) the same lognormal-type lumi nos¬ 
ity distribution with secon d-order trends dJordan et al.l 
I2007bt IHarris et al.l I2014bf) . lb) simil ar bimodal metal- 
licity distributions (jBrodie fe Straderll2006D . and (c) the 
individual GCs have King-type structures and charac¬ 
teristic scale radii that also change little with environ¬ 
ment or even galactocentric dista nce fe.g. lGeorgiev et all 
120101 1 Jordan et al.ll2005l l2007bf) . The uniformity of the 
Mccs — Mh relation can be added to this list. It appears 
to us that GC formation was controlled primarily by local 
hydrodynamic conditions in their parent GMCs and only 
secondarily by the much larger-scale and time-dependent 
structure of the protogalaxy. 

5.2. GC Metallicity 

The results from the subdivision by metallicity (Fig.Q 
are intriguing. Our main empirical finding is that GCs 
of both metallicity groups follow simple power-law forms 


Macs ~ M^ over the whole range of host galaxy halo 
mass, quite unlike the behavior of Macs versus stellar 
mass for either metallicity group. 

Considering the arguments from the theoretical mod¬ 
els described above, the strong and nearly linear correla¬ 
tion for the blue GCs is expected since the dwarf halos 
in which they formed should contain star-forming gas 
in nearly direct proportion to the halo masses, before 
feedback had severe effects (see above). But the equally 
strong correlation for the red GCs was not necessarily 
expected, and current models give less guidance for inter¬ 
pretation. The difference in slope (/3(blue) = 0.96 ±0.03 
versus /3(red) = 1.21 ± 0.03) is significant, and is the 
direct result of the trend of f re d versus Mh (see below). 

For the lower-metallicity GCs, the KG05 model gives 
a more quantitative way to understand Mccs{blue) ver¬ 
sus Mh over its entire range. In a standard picture of 
hierarchical buildup, the majority of GCs in a final large 
galaxy formed either (a) in small-to-intermediate gas-rich 
halos at early stages (z ;> 3) of the merger tree, or (b) 
in small satellites that were later accreted, but in which 
their own GCs also formed early on. In these small ha¬ 
los, the KG05 model suggests that the total GC mass 
varies nearly one-to-one with halo mass. Combining any 
number of such halos together in dry mergers without 
any additional star-forming gas will then propagate the 
same linear proportionality M G cs(blue ) ~ Mh upward 
to the giant-galaxy regime, automatically replicating the 
correlation we now see. 

The same argument does not (yet) carry over to 
the metal-richer red GCs. For these, it would not 
be valid simpl y to extrapo late the KG05 relation up¬ 
ward (as did ISpitler et all 120081) . because the metal- 
rich GCs needed to form out of more enriched gas in 
bigger halos further along t he mer ger tree . According 
to, e.g., iMuratov fc Gned3 (l2010i ): lToniml (120131) . GCs 
with metallicities [Fe/H] ~ —0.5 (the peak of the red- 
GC mode) would require host halos of masses above 
~ 10 11 Mq, several times larger than the upper limit of 
the KG05 model, in agreement with the empirical ev¬ 
idence discussed in Section 4.1 above. It would be of 
great interest to extend hydrodynamic models like KG05 
to z = 0, in which we should expect to see the rising 
proportion of metal-rich GCs as enrichment proceeds up 
to near-Solar levels in the bigger halos. 

Some metal-rich GCs may have formed in late, major 
merg ers, as can be ob s erved happening even at z ~ 0 
(e.g. I Miller et al.l 11 99 7t IWhitmore et al. 119991) . But by 
such a late stage the gas supply in the merger progen¬ 
itors is only a small fraction of the total stellar mass 
already assembled (“damp” mergers). These later merg¬ 
ers, which depend much more strongly on the contingent 
individual histories of giant galaxies, may then be re¬ 
sponsible for adding scatter at the high-mass end to the 
overall correlations in Fig. [3 

To repeat, the scaling relations (Fig. [3 and Table 
2) give M G cs{blue ) ~ M^ 96 , M G cs{red) ~ M^- 21 , 
whereas the entire population scales as Macs ~ M\ m . 
These power-law exponents have the values they do be¬ 
cause of the specific behavior of f re d or fbiue with halo 
mass, which is shown in Figure 1111 The fraction fbi ue 
exhibits a shallow power-law-like decline. Suppose then 
that we assume Mccs(blue) ~ M^ along with Mqcs ~ 
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11 12 13 14 

Galaxy Halo Mass log (M h /M 0 ) 


Figure 11. Upper panel: Correlation of the blue fraction of GCs 
versus galaxy halo mass. The solid line is the power-law fit given by 
Eq. 02 in the text. Lower panel: Red-GC fraction versus halo mass. 
The solid line is the same curve as in the upper panel, converted 
through f red — 1 - f b iue, while the dashed line is f red ~ Mj/ 19 ; 
see text for derivation. 


as before, 
we predict f Mu e 

by 


Then since f Uu e = M G cs(blue)/M G cs , 
~ = Mff 0 ' 07 . This curve, given 


fbl 


ue 


, Mh '—0,07 

1 10 10 M © J 


(3) 


and valid over the range 1O 1O M 0 < Mh <, 10 13 M©, is 
shown in Fig. [Tl] (upper panel) as the solid line. For 
Mh < 10 lo M© (dwarf galaxies) w e find em pirically that 
fbiue — 1 with rare exceptions (IGeorgiev et al.112010! ). 
while for Mh much above ~ 1O 13 M 0 the data suggest 
that the fractions may saturate at around f re d — fbiue — 
0.5. 

In the lower panel of the figure, the corresponding 
relation for f red = 1 — fbiue is shown as the solid 
line. The dashed line in the same panel shows a sim¬ 
pler power-law fit that would be expected from f re d = 
Mgcs( red)/Macs ~ M^' 21-103 = M 318 . Both curves 
provide a reasonable match within the scatter of the data. 

Lastly, we can combine the results for the red/blue 
GC subcomponents with galaxy type. In Figure |T2] we 
show the growth of the red-GC fraction f re d with galaxy 
mass, subdivided into the three approximate morpholog¬ 
ical classes. The E and SO types behave similarly. But 
interestingly, the spiral and irregular galaxies sit slightly 
above the mean line by Af red ~ 0.1, though adrnit- 
tedly the sample size is still relatively small. iKruiissen] 
(12014H predicts that galaxies experiencing fewer late 
mergers with small satellites should now have system¬ 
atically higher f re d values, since they accreted relatively 



11 12 13 14 

Galaxy Halo Mass log (M h /M e ) 

Figure 12. Correlation of the red fraction of GCs versus galaxy 
halo mass, subdivided by galaxy morphology. The solid line is the 
power-law fit to all galaxies combined, as given by Eq.[3|in the text 
and used in the previous Figure. 


fewer blue, metal-poor GCs@ In this sense the GCSs in 
S/Irr galaxies may therefore bear the imprint of a less 
active accretion history on average than the early-type 
galaxies experienced. 

In Table 2, we give the solutions further broken down 
into red and blue GCs for each of the three galaxy types, 
and these are plotted for Mccs(red) versus Mh in Figure 
M All three morphological groups separately agree well 
with the solution for all galaxies combined. For the S/Irr 
types, the higher f red fraction (Fig. [12]) largely cancels 
out the lower total Macs (Fig- EH), leaving the total 
mass in the red GCs similar to that in the E/S0 types. 

Further discussion of these points goes beyond the 
scope of this paper, but we suggest that Fig. [TT] and 
Eq. [Altogether represent one of the most important obser¬ 
vational signals from globular cluster systems in galax¬ 
ies. Where each galaxy lies on this graph is the visi¬ 
ble outcome of all the GC formation events in its indi¬ 
vidual merger tree , in all its complexity. A merger-tree 
model capable of high enough resolution to isolate GCs 
in formation and track the metallicity enrichment history 
should be able to reproduce this overall trend of the blue 
and red GC fractions. 

The intrinsic scatter present in Fig. [TT] around the 
mean relation is partly due to measurement uncertain¬ 
ties, but it should also partly be due to different individ¬ 
ual histories for galaxies of the same present-day mass. 
Two important factors in these histories are (a) the num¬ 
bers of dwarf satellites accreted at lower redshift (which 

5 More specifically, the S/Irr galaxies need to have accreted fewer 
satellites that were large enough to have globular clusters of their 
own. The graph places no limits on the accretion of very small 
satellites that may have turned into faint tidal streams but without 
any accompanying GCs. 
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lo g ( M halo/ M s) 


Figure 13. Correlation of the mass in red (metal-rich) GCs versus 
galaxy halo mass, subdivided by galaxy morphology. The dashed 
line in every panel is the fit to all galaxies combined, as given in 
Table 2. The solid line is the fit to each galaxy type individually. 


will boost fbiue)', and (b) the number of late, major merg¬ 
ers (which may boost f re d if gas is present). In addition, 
scatter may also arise because GC metallicity is not a 
strict monotonic function o f formation time. The e vi- 
dence from the Milky Way (IVandenBerg et aLll2013l) is 
that metal-richer GCs are younger by 1 — 2 Gy on the av¬ 
erage, but at any given [Fe/H] the scatter in age is at least 
1 Gy. Conversely, at a given age the scatter in metallicity 
is > 0.5 dex. The format i on models cited abo ve (KG05, 
IMuratov &; Gnediiil 120101 : iGriffen et al.1 120101) typically 
show quick buildup of metallicity in each subhalo, very 
incomplete mixing of gas, frequent merging and disrup¬ 
tion, and a rapidly changing overall potential well till 
merging abates around z ~ 2. Even as late as z ~ 3, 
the g as fractions w ithin h alos < 10 n Mg were well above 
50% (IMuratov fc Gnediiil 120101) . allowing GC formation 
in many uncorrelated starbursts over a significant metal¬ 
licity range. 


5.3. Halo and GCS Masses 

Two major criteria seem necessary to understand the 
basic correlation between Macs and Mh that we see: 

(1) Star formation is a function of very local conditions 
of gas density, temperature, and hydrodynamic state and 
is unlikely to be directly driven by the vastly larger-scale 
potential of the parent galaxy. Thus the link connecting 
the two should be the initial gas mass collected within 
the parent halos, i.e. Mqcs ~ M gas (init ) ~ Mh- 

(2) GCs during formation were relatively immune to ex¬ 
ternal feedback in their larger environments including 
UV radiation, AGN feedback, SNe, and stellar winds. 
These feedback mechanisms may hamper or truncate 
field-star formation in the large volumes of lower-density 
gas, but the GCs were forming in small regions where 


the gas density was several orders of magnitude higher 
than in their host GMCs as a whole. 

If this line of reasoning is correct, then we need to un¬ 
derstand better the structu res o f GM Cs in star-forming 
galaxies at high redshift. IKruiissenl (120141 ) has argued 
that these GMCs should be more compact than present- 
day ones at a given mass, with relatively more clumps 
of dense gas above the threshold for rapid star forma¬ 
tion. Thus in these, the GC formation efficiency (CFE, 
i.e. the ratio of GC mass to GMC mass) may be much 
higher than the present-day ratio o f — 1 % used in earlier 
d iscus sions (|Harris fc Pudritzlll994) . IMuratov fc Gnediiil 
(120101 ) also find that the fraction of stellar mass contained 
in the GCs is as high as 10% at z > 3. 

On the observational side, an increasing body of mea¬ 
surements for star-forming galaxies at z > 2 is start¬ 
ing to show (a) that molecular gas can make up as 
much as half the total gas mass; (b) that in turn, 
half the molecular gas can collect into GMCs capa¬ 
ble of quickly forming stars; and (c) the star-forming 
clouds the mselves are denser than in the present-day uni¬ 
verse (e.g. IPopping et all 120 14]: |Freundlich e t all 12014 


IDe Breuck et al. 2014 iThomson et al.l I2015D (see also 
references in lKruiissenll2014D . 

In HHH14 we proposed that the final mass ratio 77 could 
be thought of as a product of four separate ratios, 


_ / A/bary N / Mq MC \ / A/ppC N / MqC \ 

? V Mh / v Mbary ' ' MgmC ' ' AfpQC ' 

(4) 


where M bary /M h ~ 0.15 is the cosmic baryonic mass 
fraction; Mgmc/M bary is the fraction of initial gas mass 
that collects into GMCs big enough to build globular 
clusters; Mpqc/Mgmc is the mass fraction of gas in a 
typical GMC that ends up in proto-GCs (equivalent to 
the CFE); and Mgc/Mpgc is the present-day GC mass 
as a fraction of its protocluster mass. 

Of these four terms, only the first can be considered 
well known. The fourth term can be estimated roughly 
at — 0.1 given that the star-forming efficiency in a PGC 
should be ~ 0.3 — 0.5 and that the following dynamical 
evolution of the clus ter will reduce its mass by another 
facto r of 3 to 5 (e.g. lLada fc Ladal 120031 iKatz fe Ricottil 
2014). The second and third ratios, however, are espe¬ 
cially uncertain because we need to know them at high 
redshift, and because they may well show considerable 
place-to-place variance. Although far more evidence is 
needed, the works cited above are starting to point to 
the view that both {Mgmc/M bary) and {M P gc/Mgmc) 
can be of order ~ 0.1 or higher at the redshift range when 
GC formation is most active. Conservatively, we adopt 
~ 0.05 for both. 

The point of writing Eq. (4) in this form is that, even 
if some of the terms are speculative, the product of all 
four is well determined. As noted in HHH14, we can then 
invert the question: 77 can be used to constrain the second 
and third terms rather than the other way around. Our 
revised order-of-magnitude estimate, setting the product 
of the 4 terms equal to the observed value of 77 , would 
then be roughly 

77 - 0.15 x 0.05 x 0.05 x 0.1 ~4 x 10“ 5 . (5) 

Consistent with the observational work mentioned above, 
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this scaling argument indicates that several percent of all 
the initial gas mass was able to collect into GMCs mas¬ 
sive enough to build globular clusters that can survive to 
the present day. 

Lastly, of the available theoretical models that 
specifically addres s the formatio n of G Cs in hier¬ 
archi ca l merging (iDiemand et all 2005t iMoore et all 
20061: IKravtsov &: GnedinT 20051: Bekki et al.l 20081: 


Muratov fe Gnedinl 120101 : iGriffen et all 120101: iToninil 
20T1 , we find the results of KG05 to best match the data 
discussed here. For the metal-poor GCs the model suc¬ 
cessfully predicts the nearly one-to-one relation between 
Mocs(blue) and Mh , at least for halos ^3x 10 11 M Q . 
If all low-metallicity GCs form in such relatively small 
halos, then later mergers to build up a larger galaxy will 
automatically extend the correlation upward over the en¬ 
tire range of galaxy masses. The main limitations of 
the model are spatial resolution (proto-GCs are simply 
single-pixel points of high density); and redshift range 
covered (z > 3.35). Capturing the metal-rich GC popu¬ 
lation will require continuing the models to z <> 2. Mod¬ 
elling capable of going beyond both those limitations is 
increasingly possible and should open up a new range of 
tests of the theory. 


6. SUMMARY 


total initial gas mass in the halo, Alecs ~ 
M gas (init) ~ M h . 

6. An accompanying requirement for the observed 
Aloes ~ Alh correlation to make sense is then that 
external feedback processes appear to have little 
relevance for globular cluster formation. In essence, 
the relative numbers of GCs in a galaxy may then 
indicate what the global star formation efficiency 
would have looked like without the influence of 
feedback. 

7. The hierarchical model of iKravtsov fe Gnedinl 
(120051) provides a basis for understanding the ob¬ 
served correlation of the metal-poor GCs with halo 
mass, namely Mocs(blue ) ~ Mh over the entire 
range of galaxy types and masses. Extending this 
modelling to lower redshift 2 <; 3 holds out hope 
that the metal-richer GCs could eventually be fit 
into the same comprehensive picture. 

8. The specific way in which the blue (metal-poor) 
and red (metal-rich) GC fractions change system¬ 
atically with galaxy mass should provide a sensitive 
new constraint on hierarchical-merging models for 
galaxy formation at redshifts z > 2. 


We have used the database of GC system masses and 
the halo masses of their host galaxies to extend earlier 
discussions of the relation between these two quantities 
that are relics of the galaxy formation era. Our main 
conclusions are these: 

1. The observed linear correlation Moos ~ Mh con¬ 
tinues to become stronger with the growth of the 
GCS database. To first order it is insensitive to 
host galaxy morphology, applying to ellipticals, 
SO’s, spirals, and irregulars alike. A second-order 
difference between the ellipticals and spirals may 
exist in the sense that the S/Irr types were slightly 
less efficient per unit mass in forming GCs. 

2. The two GC metallicity modes show different be¬ 
haviors with Mh, but both follow extremely sim¬ 
ple scaling laws and similar low dispersion. We 
find Mocs{blue) ~ M° 96 , and (for galaxies mas¬ 
sive enough to have generated metal-ricli GCs) 
AlGCs(red) ~ M^' 21 . In short, galaxy halo mass 
provides an accurate predictor of both types of GC 
populations. 

3. The smallest dark-matter halos capable of reliably 
producing and holding metal-rich GCs ([Fe/H] > 
— 1) have Mh — 10 n MQ. 

4. We find tentative evidence that S/Irr galaxies show 
a slightly higher-than-average metal-rich GC frac¬ 
tion f re d for a given halo mass. This evidence 
may indicate that, on average, these galaxies ex¬ 
perienced fewer small-satellite accretions than did 
E and SO galaxies. 

5. We support the view that the very large-scale halo 
mass and the much smaller-scale gas-dynamical 
process of GC formation are linked through the 
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